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Abstract 

A new class of solutions of three-dimensional equations from the Boussinesq paradigm are considered. The 
corresponding profiles are not localized functions in the sense of the integrability of the square over an infinite 
domain. For the new type of solutions, the gradient and the Hessian/Laplacian are square integrable. In 
the linear limiting case, analytical expressions for the profiles of the pseudolocalized solutions are found. 
The nonlinear case is treated numerically with a special approximation of the differential operators with 
spherical symmetry that allows for automatic acknowledgement of the behavioral conditions at the origin of 
the coordinate system. The asymptotic boundary conditions stem from the l/r behavior at infinity of the 
pseudolocalized profile. A special approximation is devised that allows us to obtain the proper behavior for 
much smaller computational box. The pseudolocalized solutions are obtained for both quadratic and cubic 
nonlinearity. 
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1. Introduction: Localized Waves and Quasi-Paritcles 

John Scott-Russell observed the permanent (or "great" ) wave that can travel large distances with- 
out changing form. In 1871-1872 Boussinesq came up with a fundamental idea: the balance between the 
nonlinearity and dispersion is what makes the shape of the wave permanent [l|, 0] . He found the first ana- 
lytical solution of sech type. Later on, Korteweg and de Vries [Toj showed that in a frame moving with the 
characteristic speed, the second order in time Boussinesq equation can be reduced to a first-order in time 
equation (known now as KdV equation) for which the balance between nonlinearity and dispersion holds as 
it does for its parent: the Boussinesq equation. Boussinesq- type equations arise not only in shallow water 



flows and also in the theory of shells and plates (the famous von Karman equations) [18|, ll5| . Even the 
Schrodinger equation of wave mechanics can be shown to be a Boussinesq equation for the real or imaginary 
part of the wave function Q . Boussinesq equation is a generalized wave equation containing dispersion in 
the form of a biharmonic operator of the sought function. The dispersive effects of the biharmonic operator 
can be countered by the presence of a nonlinearity, and as a result, a permanent wave of localized type 
may exist and propagate without change. We call this balance between the nonlinearity and dispersion the 
"Boussinesq paradigm" . 



One of the main properties of the "permanent" localized waves is the collision property. Skyrme [24 
argued the idea that different particles are the localized waves of different field equations and actually found 
a particle- like behavior of the localized solution of the sme-Gordon equation [20| . Unfortunately, the latter 
solution were kinks ('hydraulic jumps') which did not appeal directly the physical intuition of particles as 
being localized 'lumps' of energy. Actually the works of Skyrme defined in essence the notion of Soliton 
which was introduced a couple of year later by Zabusky and Kruskal who discovered numerically that 
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the solitary waves of KdV equation retain their shapes after muhiple collisions. They introdueed the coinage 
soliton to emphasize this particle-like behavior. Currently the term soliton in strict sense is reserved only 
for the particle-like localized solution of fully integrable systems. For integrable and nonintegrable systems 
alike, the evolution of a wave system which is a superposition of two localized waves results in a virtual 
recovery of the shapes of the two initial localized waves but in a shifted positions: outgoing from the site 
of interaction. The collision property gave the rise of the term "quasi-particle" (or QP). When the physical 
system is described by equation(s) that conserve the energy and momentum and some integral interpreted 
as the pseudomass of a localized shape, the quasi-particles behave very much as actual particles in quantum 
physics. 

Yet there are essential differences between the quasi-particles currently known and the real particles. The 
most conspicuous is that the research is mostly concerned with ID cases, which have virtually no relevance 
to the physical reality. Even within the class of ID QPs there are more questions to be answered before they 
are shown to qualify as particles. One of the unsatisfactory features is that the quasi-particles actually pass 
through each other (in a sense, they 'percolate', in a sense, through each other) rather then scatter. This 
does not prevent the investigators to speak about scattering of solitons. To the limit of author's knowledge, 
the only work that attempts to answer the question of whether the sme-Gordon kinks percolate or scatter 



is |16l. |14||. where the Variational Approximation is applied. The standard approach without additional 



collective variables (see jl3| . and literature cited therein) cannot answer this question. 

It has been just recently discovered (see [2^ and the further corroboration in [2^) that under certain 
conditions the solitons of the System of Coupled Nonlinear Schrodinger Equations (SCNLSE) (called alter- 
natively Vector NLSE) can actually scatter without crossing each other (percolating through each other). 



1.1. Classical Solitons as Quasi-Particles: The Con's 

Mathematicians place emphasis on the full integrability of the system that exhibits soliton solutions. 
However, physically speaking, the three laws known to govern motion are the conservation of mass, energy, 
and momentum. The full integrability is not directly relevant for the physical implication of localized 
solutions as quasi-particles. 

There are essential differences between the real particles and ID quasi-particles known currently for 
which hundreds of strict mathematical results are proven. First, the known quasi-particles actually pass 
through each other (percolate in a sense), i.e. they do not repel (scatter) like the real particles. Second, 
the repelling or attracting of the QPs depend very much on the potential of interaction defined by the 
asymptotic behavior of their tails. In ID, for virtually all main soliton supporting equations the potential 
is attractive, and decays exponentially with the distance between the QPs. In real-world physics, at large 
distances the potential (producing the gravitation force) decays algebraically as the inverse of the distance, 
which is a radical difference from an exponential decay. Third, the results about the interaction in more 
than one dimension are rare and mostly qualitative, not quantitative. The repelling or attracting of the QPs 
depend very much on the potential of interaction defined by the asymptotic behavior of their tails. 

The above list of discrepancies between the currently known QPs and the actual physical particles can 
be extended, but the main points are those mentioned above. 



1.2. The First Forays in the Realm of 2D Boussinesq Solitons 

The first thing that comes to mind is that this exponentially decaying asymptotic behavior is due to 
the one-dimensional nature of known QPs. Unfortunately, very little is known for solitary waves that are 
truly localized in the two spatial variable. For the time being, the main success has been in finding solitary 
waves of the Kadomtsev-Petviashvily (KP) equation (see the original article 17 1, also [12, HH l22l | and the 
literature cited therein). 

The strictly localized steadily propagating solutions of the truly 2D equations of Boussinesq type were not 
investigated until recently. The question is very important and we have developed three different numerical 
techniques to interrogate the problem. A very fast Galerkin spectral method is used in Q ; a semi-analytical 
method for relatively small propagation speed is presented in [lO| ; and a finite-difference solution with special 
implementation of the asymptotic boundary condition is given in Q . The three methods agree quantitatively 
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very well in the common domains of the values of the parameters, so that the findings can be considered 
validated with high accuracy. 

The findings about the 2D localized solution of Boussinesq equation can be summarized as follows. The 
profile of the standing soliton has a super exponential decay at infinity, namely e""^ /y/r. This property is not 
robust, and even for very small propagation speeds, the decay at infinity changes to algebraic, proportional to 
However, the most disappointing property is the fact that the steadily propagating 2D shapes are not 
stable and after some time they either dissipate or blow up. This important result was first found by means 
of a special new scheme Q and confirmed quantitatively in [ll| with a different numerical approximation. 
All these findings rule out the solution from @, S 0, Eli ^ ^ candidate for the model of the real particles. 
The instability in time may be tackled by adding specially selected nonlinear terms that can either alleviate 
the blow-up or prevent the dissipation. The non-robustness with respect to the parameter giving the phase 
speed of the structure tells us that in order to make the standing localized wave moving, the profile behavior 
up to infinity has to be instantaneously changed. 

1.3. The Notion of a Pseudolocalized Solution 

The disappointment with the, so to say, "strictly" localized solution requires that we reconsider the 
notion of localization. Since the equations that describe the field for a given family of particles are still 
under construction, one cannot be a-priori sure what is the meaning of localization and what is the meaning 
of the primary variable. The requirement for strict localization that the profile function should be square 
intcgrable, while a closer look reveals that the quasi-particle behavior may be well satisfied if merely the 
Laplacian (the curvature) or the gradient of the profile (the slope) are square intcgrable. Such kind of 
solutions were first considered in [7[ and called "pseudolocalized" . 

The idea of the pseudolocalized solution is that one solves a nonlinear equation for the strictly localized 
curvature and then integrate the Poisson equation with r.h.s being the known curvature to obtain the profile 
of the original amplitude. In 2D, such an integration will create a logarithmic behavior at infinity. Fortu- 
nately, in 3D the behavior of the solution of the Poisson equation is as 1/r, which is very agreeable physically. 
Despite of the visual resemblance to a strictly localized solution, the profile is actually pseudolocalized in 
the sense of the above definition. 

The ultimate question to answer is whether the 3D solitons can retain their shape when they are allowed 
to evolve according to time dependent equation. To demonstrate the presence of the collision property 
numerically in 3D seems like an insurmountable task due to the required computational resources. The 
feasible goal is to first show that the shape of the soliton does not change when left to propagate according 
to the full equation. A prerequisite to investigating the collision properties of the pseudolocalized solutions 
is to find analytical or numerical representation of a single QP. This is the goal of the present paper. 

2. Selecting the model 

As already mentioned, we focus on the Boussinesq equation which arises in shallow water flows and 
in the theory of shells and plates. It is a generalized wave equation containing dispersion in the form of 
a biharmonic operator of the sought function. The dispersive effects of the biharmonic operator can be 
countered by the presence of a nonlinearity, and as a result a permanent wave of localized type my exist 
and propagate without change 

Consider a generic version of Boussinesq equation 



Alternatively, in shell theory the nonlinearity enters the picture in a different fashion (see the simplified 
version from [tI) 



utt = A[u ~ oi2i^ — OLz"^ ~ /3Au]. 



(1) 




(2) 



X = p cos B cos 



y = p cos 9 sin 



z = psm 
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and restrict ourselves to the ease of spherieal symmetry when 

A - — — ^ — 

dp^ dp 

3. The Linearized Equation and the notion of Pseudolocalized Solution 

The hnearized version of the equation is merely the Euler-BernouUi equation, i.e., 

utt = Aw - /3A^u. (4) 

In 2D, the solution with point symmetry that does not exhibit a singularity at the origin was first found 
in different eontext (see, e.g., Q) to be 



u 



(2)=Ao(r/v/^)+ln(r/v/^), (5) 



where Kq is the modified Bessel function of the second kind. Clearly, this solution is not localized function, 
because it increases to infinity at r oo. What is is more important is that it is not even pseudolocalized, 
in the sense of the definition of the present work, because the following integral diverges (denote 77 = -^)'- 



/ [uf'^frdr=l [{K'Q{T])f ri + 2K'Q{ri) + const -hYe + 111 R, for i?>l. 



(6) 



Note that the integral diverges also at e — > 0. 

In 3D, the solution of the linear problem with point symmetry that does not have a singularity at the 
origin is (see 0]) 

,(3)(^)^l-exp(-p/V?)^ 
P 

which decays at infinity, but is not localized, because the integral of the square of it diverges linearly at 
infinity. However, it is a pseudolocalized solution because the L^-integral of the gradient converges, namely 

M^Kp) .2 2, _r\ l-cxp(-p/V?) exp(-p/V?) i2 _ 1 

T P '^P — 2 ' 75 P '^P — TTTS- \°) 

op Jo L p2 J 2y/P 

It is important to mention here that the above solution does not exhibit a singularity at the origin, but 
it is not a smooth function of the coordinates (x, y, z). It appears to be akin to the "peakons" from [3| (see, 
also the in-dept study of the numerical approach to peakons evolution Q )with the difference being that it 
decays at infinity as 1/p rather than e"''. In addition, we have not been able to show that an analytical 
expression for the peakon can be obtained when the nonlinearity is not neglected. Actually, this paper deals 
with finding a numerical solution for the nonlinear case. An important advantage of the pseudolocalized 
solution is that if it propagates, then the additions to the profile is of order I/r^, which does not qualitatively 
change the decay at infinity, as it is the case of strictly localized solutions. The profile stays in the the same 
class of localization. 



4. Effect of Nonlinearity on the Shape of the Pseudolocalized Solution 

The most conspicuous trait of the above defined pseudolocalized solution is that it is not smooth at the 
origin of the coordinate system. It is interesting to note that for some specific nonlinearity, such as the 
one in Eq. the effect of nonlinearity leads to smoothing of the profile at the origin (see 0), and the 
pseudolocalized solution does not have the peak of the linear approximation. However, the way nonlinearity 
appears in Eq. ^ does not act to smooth the peak. There is no known analytical solution for the nonlinear 
version of the governing equation, so we present here a numerical solution. 
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The full three-dimensional model requires significant computational resources. For this reason we have 
already restricted ourselves to the case of spherical symmetry. Now, the sought function depends only on 
the radial coordinate p. To emphasize this fact, we use a new variable /(p) := u{x,y,z). Then the original 
equation is a equivalent to the following system 

/3^^P^^/ + «2/^ + «3/-^-/ = 5, (9) 

(10) 

dp dp 

One can integrate the second of the above equations to obtain g = Ci/ p + €2- The truUy localized 
solution requires that Ci = C2 = 0. The pseudolocalized one is obtained when only C2 = 0, and Ci can 
have an arbitrary nontrivial value. We use this fact to reduce the problem to solving numerically the single 
equation 

P\^P'^f + »^f + ^3f~f=-^ (11) 

dp dp p 

where a is at our disposal and choosing it will lead to different solutions. The boundary conditions are the 
behavioral condition at the origin, and the asymptotic condition at p = Poo which requires that 



dp Poo 

where poo is the finite value of the radial coordinate at which the computational domain is truncated. 



4-1- Difference Scheme 

The boundary value problem Eqs. ([TT|) . ((T^ is to be solved numerically. We use a grid which is staggered 
by from the origin p — Q. while it coincides with the "numerical infinity", p — poo. Thus 

Pi^^ii-D^^ Pi±^=Pi^\h, i = l,2,...N, 

where h — Poo/{N — 0.5) and is the total number of points. 

The staggered grid gives a unique opportunity to create difference approximations for the Bessel operators 
involved in our model that take care of the singularities of the respective Bessel operator automatically, 
without the need to impose explicit behavioral boundary conditions in the origin. This is made possible by 

the fact that 01 = 0. 

2 

The main difference of the above and the smooth (non-peakon) solution is that the latter is a result 
of a bifurcation and exists only for one specific value of the amplitude of the soliton. The pseudolocalized 
solution exists for continuous interval of the parameter a in Eq. 

The finite-difference complete scheme with the boundary conditions is given by the following system with 
tree-diagonal matrix: 

4/2 - [p| + h'pflfi = -h'pl[a2{f^)' + «3(/r)'] + ah^p,, (13a) 
P-+J^+i - k'_i + + h'pj]f^ + plJr-i = -h'pl[a2{f:f + a,{f:f] + ah^p., 

i = 2,...,N-l, (13b) 

In+i -{— + l)fN = 0. (13c) 
Pn 

Note that the first line is the so-called behavioral conditions stemming from the properties of the operator 
rather than imposting additional boundary condition. 

To control the convergence we use also relaxation /"^^ := a;/ + (1 — w)/", and repeat the iterations 
until convergence is attained when the value of the last iteration is taken as the solution for function /. 
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4-. 2. Results and Discussion 

All numerical experiments to follow are carried out with P ~ 1, since one can use a simple rescaling of 
the independent spatial variable. We set poo = 1000 and use N ~ 20001 for the total number of grid points. 
In most of the cases reported below, a relaxation parameter uj — 0.2 was sufficient to secure convergence of 
the iterations. 

First we examine the case of only quadratic nonlinearity, i.e., 012 = 1 and a^, = 0. The interesting 
observation is that the algorithm is not convergent for negative a when the profile becomes positive. For 
positive a, the profile becomes taller with the increase of a (see the rightmost panel of Fig. [T]), but the 
increase is lesser than the increase of the coefficient a as seen from the left panel of Fig. [1] where the profile 
is scaled by a. The middle panel shows that for very large distance from the origin, all of the scaled profiles 
collapse on the linear solution. One can say that the nonlinearity acts to reduce the overall amplitude of 




2 4p 6 8 10 1 10 p 100 1000 2 4p 6 8 10 



Figure 1: The peakon solution for tlie quadratic nonlinearity Q2 = 1 and 03 = and for different values of parameter a from 
Eq. mill . In the left and middle panel (log- log plot), the profile is scaled by a. In the right panel is the original profile. 
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Figure 2: The peakon solution for the cubic nonlinearity 02 = and as = 1 and for different values of parameter a from 
Eq. I lllll . In the left and middle panel, the profile is scaled by a. In the right panel is the original profile. 



the solution, all other conditions equal. 

In a similar fashion we compute the numerical solution for the cubic case 02 = and a2 7^ 0. It is well 
known that for the ID cubic Boussinesq equation, the case > does not lead to localized solutions. More 
specifically, for > solitons of sec/i-type cannot exist. It is interesting to note that we actually found 
pseudolocalized solutions even for this case, but only for a < 0.914. In fact, for the last case the solution 
was obtained only when we reduced the relaxation parameter to w = 0.005 and increased the number of 
grid points to 40001, all other parameters being the same. However, all our efforts to obtain a result for 
a = 0.915 did not meet with a success. We also mention here that for cubic nonlinearity, the positive and 



6 



negative amplitudes are equally possible for the solution, and we present here the case of a > which results 
in negative amplitude of the peakon. The results are presented in Fig. [51 Now, as it should have been 
expected, the profile is lower than the linear one for all values of a (see the right panel of Fig. (2). Yet, the 
nonlinearity retains its steepening effect (see the left panel of Fig. [2l). Naturally, because of the smaller 
values of a, the steepening effect is less obvious. 

Finally, we treat the cubic case with as < 0, when the energy functional is not positive definite allowing 
the bifurcation and appearance of ID sech solutions in the ID case of strictly localized solutions. We found 
solutions with very large a. With lo = 0.2 we easily reached a = 18, but for larger a the relaxation parameter 
must be smaller. The result suggests that eventually, the solution may exist for arbitrary a, but to claim 
this, one needs an analytical result. Once again, both positive and negative solutions are possible and are 
equivalent upon changing the sign of the amplitude (or what is the same: the sign of a). Fig. [3] presents 
the result for negative a, when the solution has a positive amplitude. Qualitatively, it is very similar to the 
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Figure 3: The peakon solution for the cubic nonhnearity 02 = and 03 = —1 and for different values of parameter a from 
Eq. Illlll . In the left and middle panel, the profile is scaled by a. In the right panel is the original profile. 



quadratic case, with the small quantitative difference that the profiles are now steeper for the same a than 
in the quadratic case. 



5. Conclusions 

In the present work we have argued the need to look beyond the strictly localized solutions when the 
concept of Quasi-Particle (QP) is applied. We consider the Boussinesq equation as the featuring model, 
and introduce the concept of pseudolocalized solution which does not belong to the L^-space itself, but its 
gradient and Laplacian do. It turns out that the pseudolocalized solution exists only in dimension three or 
greater. A similar attempt in 2D produces a clearly no-localized profile which behaves at infinity as Inr, 
where r is the polar coordinate. Its gradient behaves as l/r, which means that it is not localized either. 
Thus the 2D solution is not pseudolocalized. 

For the 3D case we found a pseudolocalized solution for the linearized version of the equation and showed 
that its profile decays at infinity as 1/ p (where p is the spherical radial coordinate), which is rather different 
from the exponential decay of the strictly localized solution. We shown that the profile is not smooth at the 
origin of the coordinate system, and can be called "radial peakon" . 

We devised a difference scheme to solve the full nonlinear equation numerically and showed that retaining 
the nonlinearity does not change qualitatively the profile obtained from the linear equation. Naturally, at 
infinity the profile decays with the same law 1/ p, since the nonlinear solution gradually approaches the linear 
one. The nonlinearity plays a role in steepening the peak of the solution at the origin. 

The importance of this new kind of solution brings the concept of quasi-particle a step closer to modeling 
real particles in physics as localized solutions of a field equation (see the pioneering works of Skyrme [23.[20l|V 
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since the potential of the interaction of two QPs is defined by the decay of their "tails" . Since the decay 
here is the predominant part of the asymptotic behavior of the potential of interaction will be close 
to 1/p which yields to a gravitation-like force of interaction between the radial peakons. The results of the 
present work may prove instrumental in a new line of investigation in which the QPs are represented by the 
radial peakons investigated here. 
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